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Abstract. We extend a coherent network data-analysis strategy developed earlier 
for detecting Newtonian waveforms to the case of post-Newtonian (PN) waveforms. 
Since the PN waveform depends on the individual masses of the inspiraling binary, the 
parameter-space dimension increases by 1 from that of the Newtonian case. We obtain 
the number of templates and estimate the computational costs for PN waveforms: 
For a lower mass limit of IMq, for LIGO-I noise, and with 3% maximum mismatch, 
the online computational speed requirement for single detector is a few Gflops; for a 
two-detector network it is hundreds of Gflops and for a three-detector network it is 
tens of Tflops. Apart from idealistic networks, we obtain results for realistic networks 
comprising of LIGO and VIRGO. Finally, we compare costs incurred in a coincidence 
detection strategy with those incurred in the coherent strategy detailed above. 
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1. Introduction 

Close compact binaries are among the prime sources of gravitational waves that hold 
promise for detection with upcoming laser interferometric detectors such as LIGO, 
VIRGO, GEO-600, TAMA, and AIGO. The back-reaction of radiated gravitational 
waves results in an inspiral with an eventual merger of the two companions of the binary 
system. This adiabatic inspiral waveform has been accurately modeled upto 2.5 post- 
Newtonian order (PN) In an earlier work [^, we developed a formalism for detecting 
inspiral waveform with a network of detectors. The proposed analysis is of coherent 
nature where the network is treated as a single detector and the data is combined using 
the phase information optimally. In 0, we used the maximum likelihood detection 
(MLD) technique, which involves correlating the output of a network of detectors 
with the family of expected waveforms (or templates) and selecting the maximum of 
the network likelihood ratio for decision making 0. To reduce computational costs 
involved in searching over the space of source parameters, we succeeded in analytically 
maximizing over 4 of these parameters, namely, the overall amplitude, initial phase and 
the orientation angles of the binary orbit. The maximization over the time of arrival 
(or, alternatively, over the time of final coalescence) of signal was carried out via FFTs. 
Estimates of computational costs involved in searching over the source-direction angles 
and the chirp mass were obtained for the simplistic case of Newtonian waveforms. In 
this work, we extend the coherent network analysis to the more realistic case of PN 
waveforms. A restricted PN waveform depends on individual masses of the companions 
instead of the combined chirp mass. This increases the number of parameters by one. 
We estimate, in general, the costs involved in searching over the masses as well as 
the source-direction angles for realistic network configurations. Finally, we describe a 
coincidence network detection strategy and compare costs incurred in it with those in 
the coherent detection strategy. 

2. Restricted post-Newtonian signal at the network 

The signal s\t) at the constituent J-th detector of the network is given by 



where k is the overall amplitude that depends on the fiducial frequency fs and the 
masses of the binaries. 6c is the phase of the waveform at the time of final coalescence. 
The extended beam-pattern functions of the J-th detector, Ej, depend on the source- 
direction angles, {9, 0}, the orbit orientation angles, {e, ip}, the J-th detector orientation, 
(«(/),/?(/), 7(7) ), and the sensitivity, g(^j), of the detector to the incoming signal. Finally, 
S\t) is a normalized complex signal such that in the stationary-phase approximation 
(SPA) its Fourier transform (FT) is 
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for positive frequencies. Above, the phase of the 2.5 restricted PN waveform at the J-th 
detector is the scalar 

vI/(,)(/;/„t„M,r/,n3,ni) =t9^e(/)M(/;/s), (3) 

with the parameters -(9^ consisting of the final-coalescence time tc, the total mass M, 
the mass ratio t]{:= 11111712/ M'^) and the source-direction described by two components 
ui and ns of unit vector n pointing to the source. Given below are and ^(i)^ '■ 
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Here, and are, respectively, the x and z coordinate values (in units of c/ fg) of 
the location of the J-th detector in a fiducial reference frame. 



3. Number of templates 

In this section, we estimate the number of templates required to search over the 
parameter space. In Ref. we showed that for a given pair of source-direction angles 
{0,(f)), the network likelihood ratio, when maximized over the overall amplitude, Sc, e, 
and ip, gives the network detection statistic. Numerical maximization of the statistic 
over the rest of the parameters, namely, masses and source- direct ion angles is performed 
by using a template bank. We estimate the number of templates by calculating the 
volume of the parameter space of interest obtained by the computing metric on the 
manifold and dividing by the size of each template. When the network statistic is 
dependent on the parameters solely through the difference between the parameter values 
of the signal and the template, then the metric on the parameter manifold is flat and, 
hence, the template placement is uniform. 

It is well known that with PN order > 1, the metric on the manifold is not fiat. 
The Tanaka-Tagoshi ||^ coordinates provide a convenient and an elegant way to carry 
out further analysis. The salient feature of these coordinates is to make the metric 
Euclidean on a fiat manifold, which is an approximation to the actual manifold. Also 
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the coordinate volume of the parameter space in these coordinates is same as the proper 
volume which immediately gives the number of templates as has been described above. 

4. Computational costs 

In this section, we estimate the cost involved in numerically searching over the rest of 
the parameter space mentioned above. This cost has two important components: 

(i) The cost involved in FTs : MLD technique requires to cross-correlate the data with 
all possible templates in the rest of the parameter space involving mass parameters 
and the direction angles. Since information of the direction angles is encoded in 
time delays, network correlation vectors for templates differing in direction angles 
can be constructed by combining the correlation outputs from different detectors 
with appropriate time-delays as described in 0. Thus, the cost involved in FTs 
is equal to the number of computational operations required in searching over the 
intrinsic parameters, in our case, the two masses of the binary. 

(ii) The cost involved in scanning the time-delay window: The optimal statistic needs 
to be evaluated by combining the correlation vectors with appropriate time-delays. 

Consider a network of Nd detectors. Let be the number of sampled points in 
a data train at each constituent detector. If the templates are stored in memory, then 
the computing cost in FTs is 6NDNnxi-x2 log2 where nxi-x2 reflects the number of 
templates in mass parameters. The number of floating point operations to construct a 
network statistic for one pair of direction angles is SN^,. If Utot is the total number of 
templates then the total computational cost is 

C = 2NDN{8ntot + ^nx.-x^ loga N) . (5) 

Online data processing requires the data processing rate should be equal to the data 
acquisition rate. Thus the length of the data which is effectively processed is equal to 
the length of the zero padding. We obtain the online computational speed by dividing 
the cost by the length of the padding interval. We use the analytical fits to the noise 
curves of LIGO and VIRGO which we enlist in Table We tabulate the results for 

Table 1. Analytical fits (for positive frequencies) to noise power spectral densities, 
2sh{f), of the interferometric detectors studied in this paper [Q. We take Sh{f) to be 
infinite below the seismic cut-off frequency fs- We choose the high frequency cut-off, 



to be 800 Hz for aU /. 


Detector 


Fit to noise PSD, 10^^ x s/,(/)/Hz-^ 


fo (Hz) 


fs (Hz) 


VIRGO 


3.24 [(6.23///o)-^ + 2(/o//) + 1 + {f / fo?] 


500 


20 


LIGO I 


9.0 [(4.49///o)-56 + o.l6(///o)-4-52 + 0.52 + 0.32(///o)2] 


150 


30 



various idealistic as well as realistic networks in Table |^. For the case of real network 
of LIGO- VIRGO P] with their respective noises, we estimate the average number of 
templates for most of the astrophysical range of e and ip to be Utot ~ few times x 10^° 
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Table 2. Number of templates, computational costs, and online computing speeds 

required for a search using specific networks. The detector networks are labeled as 
I for a single detector. III for three identical detectors with identical orientations 
placed on Earth's equator forming an equilateral triangle. The detector Xb denotes 
a detector with LIGO-I noise at the location of the detector D. The letters L, H, V, 
T, and A denote, LIGO detector at Louisiana, LIGO detector at Hanford (of 4 km 
arm-length), VIRGO, TAMA and AIGO sites, respectively. We assume LIGO-I noise 
for both the LIGO detectors. We present results for lower mass limits of 0.5Mq and 
I.OMq. The maximum length of the 2.5 PN chirp is 96.5 sees and 306 sees for minimal 
mass limits of IM0 and O.BMq respectively. We assume fiducial frequency /, = 30 
Hz except for the LV case. We consider data trains of 1100 sec. for O.5M0 and 400 
sec. for I.OM0 sampled at 2 kHz so that N ~ 10^. For the LV network, the length of 
the longest chirp is ^ 284 sec. for I.OMq and ^ 900 sec. for O.5M0. The number of 
points in the data train ~ 10^ — 10''. The mismatcli is taken to be 3%. 



Network 


mass limit 


ntot 


nxi-X2 




S 


configuration 


(Me) 


(xlO^) 


(xlO^) 


(fl-pt ops) 


(Gflops) 


/ 


0.5 


0.0214 


10.7 


0.3 


0.37 




1.0 


0.0042 


2.1 


0.02 


0.06 


LH 


0.5 


0.81 


16 


37.5 


4.7 




1.0 


0.16 


3.2 


0.37 


1.2 


LXv 


0.5 


2.2 


16 


8.7 


19 




1.0 


0.44 


3.2 


0.72 


2.35 


LXx 


0.5 


2.67 


16 


10 


12 




1.0 


0.52 


3.2 


0.82 


2.7 


LXa 


0.5 


1.7 


16 


6.8 


15 




1.0 


0.68 


3.2 


1 


3.36 


LV 


0.5 


14.1 


93 


160 


67 




1.0 


2.7 


19 


11 


12 


III 


0.5 


1.5 X 10^ 


16 


7.7 X 10''^ 


9700 




1.0 


2.9 X 10^ 


3.2 


5.5 X 10^ 


180 



for lower mass limit of 0.5Mq. We take data trains of length 3000 sec. corresponding to 
the longest chirp of ~ 900 sec. for VIRGO. Taking a samphng rate of 2 kHz., the data 
must be processed in 2100 sec. The online data processing demands a computational 
speed of few thousands Tfiops. 

Here, we note that for networks with No > 3, the computational cost required 
to construct optimal network statistic while searching over the source-direction angles 
overshoots the FT costs. As a result, the computational requirements are beyond the 
reach of the current technology for a flat search. 

5. Coincident Search 

To focus on the essential aspects of a coincidence search strategy, we consider 
the simplistic case of a network comprising of detectors with identical noises and 
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orientations, but with arbitrary locations. The network detection statistic in such 
a search is taken to be the minimum element in {|C^|, |C^°|}, where is the 
single detector statistic evaluated on the data of the J-th detector (see Ref. 0). 
Therefore, unlike a coherent search, a coincident search involves first establishing 
threshold-crossing by the single-detector detection statistic in each of the detectors 
in a network. Furthermore, claiming a detection by the network requires that the 
parameters corresponding to the threshold-crossing templates lie within error intervals 
of one another, such that they can be consistently ascribed to a a single astrophysical 
event. This requirement alone immediately implies that, even in this simplistic network 
the computational cost in a coincident search is larger than N^Ci, where Ci is the 
computational cost for a single detector search. 

To ascertain exactly how large this cost is, we first describe for a network of two 
identical detectors our search algorithm, which is based upon a most powerful search 
and is not necessarily the cheapest computationally 0: 

(i) Filter the data x^(t) and x^(t) from the two detectors, respectively, with a bank of 
single-detector templates to draw two separate lists of threshold crossers. Label 
these "candidate events" = E{tf]'d]) and = E{tf]^]), respectively, 
where tf^ (t^^) denotes the time of arrival of event i at detector 1 (2), and 
i = l,2,...,m, j = l,2,...,n. Note that m 7^ n, in general. Also, -d] denotes 
the template-parameter vector characterizing event i at detector 1. The above 
nomenclature is suited to handle the possibility of two or more templates triggering 
off simultaneously, say, on the data from detector 1. In such a case, one will have 
more than one event with t^l^ = t^^ = t^l^^, but with 7^ "d] 7^ 

(ii) "Time window" veto: Let detector 1 have a smaller number of candidate events 
than detector 2. With each E} associate a set W'^(t^^; T^\j) of candidate events Ej, 
such that f} — Tl\j ^ tj < t°; + Tl\y Here, tI^^ is the sum of the light-travel time 
between the two detectors and the sum of magnitudes of the estimated errors in 
their arrival times at detectors 1 and 2. Note that an event E'j may appear in more 
than one set. That is, it may happen that E'j G Tl\j)nW'^{t'^] '^1%)-, where 
i ^ k. Discard from the lists those E"^ that do not belong to any W'^itf^; tI\-). 

(iii) "Parameter window" veto: Compute the covariance matrix in the parameter space 
around E} and around each event in W'^itf^] ^1%) from the ambiguity function ||^. 
Estimate the parameter error, (A'd^j), to be the square-root of the variance 



of the parameter '&] (i?^) derived from this matrix. Discard events in W'^{t1^; 
that have ["i? • — ■(9^ '^l > jA-i?- '^I + {A-t^^ ^\ for each parameter index /i. 



12 ^ 

c ij, 



The pairs of candidate events surviving the above vetoes are the "detected" events. A 
more sophisticated approach involving further vetoes of the type discussed in Ref. 
will be studied elsewhere. 

The above steps explicitize the computational costs, over and above that of N^Ci, 
that are necessary in a coincidence detection, but are often glossed over: Extra costs are 
involved in computing parameter errors and implementing vetoes based on them. These 
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costs obviously scale as the number of the candidate events in each detector (whereas, 
the cost in a coherent search is independent of it). These counts, in turn, depend on 
the value of the detection threshold and, therefore, on the false-alarm probability. The 
number of floating point operations (Flop) needed to estimate the error in a parameter, 
^, is close to that involved in taking the second derivative of with respect to 
^ . Using the discrete version of second derivative the number of Flop involved 
~ lOCi. Therefore, in an 8 dimensional parameter space (based on the independent 
parameters (r, (5c,'(9°, ...,^9^), the number of Flop required to estimate parameter errors 
for all candidate events is about 80Ci x Y^^=i ^ij where Nj is the number of candidate 
events in detector I. Additional operations required to compare the parameter values 
across detectors (using the inequality given in step (iii) above) and veto events scales 
as n7=i Nj, which is a small fraction of the total cost for Nd = 0{1) and Nj = 0(10). 
Thus, neglecting this last contribution, the total number of Flop scales as: 

Nd 

NdCi + 80Ci Y.^! ■ (6) 



7=1 

l2 ;„ /vr ^ in6 



For comparison with the coherent search costs, we take Nj = 10 in ~ 10 data 
points in each of the 3 detectors in a network. For, a minimum mass of O.SMq, Table 2 
shows that Ci = 0.3 x 10^^. Thus, for network configuration III, the total number of 
Flop in a coincident search is about 7.2 x 10^^ which is very close to Ctot = 7.7 x 10^'' for 
a coherent search. One may argue that it is possible to reduce Nj in each detector by 
using additional vetoes of the type adopted in Ref. flOl- Such steps will surely reduce 



the contribution from the second term in eq.(P). Nevertheless, the additional costs in 
implementing such vetoes is very large as well and must be explored in more detail. 

It is easy to see that with more events or more detectors, the cost related to Eq. (P) 
can only rise. With future detectors, where the detection thresholds will be lower than 
present ones owing to their higher sensitivity, the number of events with larger signal- 
to-noise ratios will increase, consequently, increasing the computational cost further of 
a coincidence search. 



6. Conclusion 



As shown in Table 2., the computational cost in a coherent search rises markedly in 
going from Nd = 1 to 3. This is expected because the number of parameters and, 
therefore, the parameter volume accessible to a search increases from 5 (for Nd = 1) to 
9 (for Nd = 3). Indeed, for Nd > 3 the cost required to search over source-direction 
angles overshoots that required for the FFTs. Beyond Nd = 3, the computational cost 
in a coherent search, however, stabilizes. This must be contrasted with the cost behavior 
in a coincidence search, where it continues to increase with Nd- Specifically, given a 
network of identical detectors and a false-alarm probability, for a low enough detection 
threshold, a coincidence search will cost more than a coherent search for Nd > 1. In 
either search, the computational costs are very large and, hence, call for investing in the 
exploration of more efficient search techniques, such as hierarchical strategies. 
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